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ABSTRACT 

We examine the effects of a relativistic wind on the spin down of a neutron 
star and apply our results to the study of Soft Gamma Repeaters (SGRs), 
thought to be neutron stars with magnetic fields > 10^^ G. We derive a 
spin-down formula that includes torques from both dipole radiation and episodic 
or continuous particle winds. We find that if SGR1806-20 puts out a continuous 
particle wind of 10^"^ erg s~^, then the pulsar age is consistent with that of 
the surrounding supernova remnant, but the derived surface dipole magnetic 
field is only 3 x 10^^ G, in the range of normal radio pulsars. If instead, the 
particle wind flows are episodic with small duty cycle, then the observed period 
derivatives imply magnetar-strength fields, while still allowing characteristic 
ages within a factor of two of the estimated supernova remnant age. Close 
monitoring of the periods of SGRs will be able to establish or place limits on 
the wind duty cycle and thus the magnetic field and age of the neutron star. 
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1. INTRODUCTION 



There is a growing collection of pulsating high-energy sources which occupy a unique 
phase-space in their combination of long (> 5 s), monotonically increasing periods and 
high period derivatives. One subgroup in this collection are the soft gamma-ray repeaters 
(SGRs), transient sources that exhibit repeated bursts of soft (~ 30 keV), short duration 
(0.1 s) 7-rays. Bursts of average energy 10^^ ergs repeat on irregular intervals, while giant 
bursts of energy ~ 10^^ ergs have been observed once in each of the sources SGR 0526-66 
and SGR 1900+14. Recently, a period of P = 7.47 s has been detected from SGR 1806-20 in 
quiescent emission (Kouveliotou et al. 1998) and a period of P = 5.16 s from SGR 1900+14 
in both quiescent (Kouveliotou et al. 1999) and giant burst emission (Hurley et al. 1999a). 
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Both have measured period derivatives around P ~ 10^^'^ss^^. Another group of sources 
having similar P and P are the anomalous X-ray pulsars (AXPs), pulsating X-ray sources 
with periods in the range 6 — 12 s and period derivatives in the range 10~^^ — lO^^^ss"^ 
(Gotthelf & Vasisht 1998). These sources have shown only quiescent emission with no 
bursting behavior. 

The most plausible and commonly invoked model to explain the characteristics of both 
the SGRs and AXPs is that of a neutron star having a dipole magnetic field of > 10^"^ 
G, much higher than the fields of ordinary, isolated pulsars and well above the quantum 
critical field of 4.4 x 10^^ G. Such stars, known as magnetars, were first proposed by Duncan 
&; Thompson (1992), Usov (1992) and Paczynski (1992) to account for various properties 
of SGRs and 7-ray bursts. All of the SGR and AXP sources are found in or near young 
(r < 10^ yr) supernova remnants (SNR). SGRs in addition have been associated with 
X-ray and radio plerions whose emission power far exceeds the dipole spin-down power. 
It was therefore proposed that relativistic particle outflows from the SGR bursts (Tavani 
1994, Harding 1995, Frail et al. 1997) or from a steady flux of Alfven waves (Thompson & 
Duncan 1996), provide the power for the nebular emission. The existence of such a wind 
has been inferred indirectly by X-ray and radio observations of the synchrotron nebula 
GlO.0-0.3 around SGR1806-20 (Murakami et al. 1994, Kulkarni et al. 1994)[[ Thompson 
& Duncan (1996) estimated that the particle luminosity from SGR1806-20 is of the order 
of 10^^ erg s~^. Such energetic particle winds will also affect the spin-down torque of the 
star by distorting the dipole field structure near the light cylinder (Thompson & Blaes 1998 
[TB98]). We will show, however, deriving a formula similar to the one given by TB98, that 
if relativistic wind outflow continuously dominates the spin-down of the neutron star in 
SGR 1806-20 at a level ~ lO^^ergs"^, then the surface magnetic dipole field is too low to be 
consistent with a magnetar model. It is possible however that the wind outflows from SGR 
sources are episodic, lasting for a time following each burst that may be small compared to 
the time between bursts. In this case, dipole radiation will dominate the spin-down between 
bursts, even though wind outflow may dominate the average rotational energy-loss rate. 
Using a general formula for the spin-down torque that includes both dipole radiation and 
episodic particle winds, we derive the magnetic field and characteristic age of the neutron 
star as a function of the observed P and P, the wind luminosity and wind duty cycle. We 
find that the derived surface field and age can have a range of values between the pure 
dipole and pure wind cases that depend on the duty cycle of the wind outflows, even for a 
constant value of the average wind luminosity. 



^Note that a revised IPN location now places the SGR source outside the core of the radio plerion (Hurley 
et al. 1999b) 
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2. ENERGY LOSS FROM A WIND 

It is known that neutron star rotation drastically distorts the magnetic field near and 
beyond the light cylinder radius i?LC = c/VL (where f2 = 2ti/P is the neutron star angular 
velocity; P is the period), and that magnetic field lines crossing the light cylinder remain 
open all the way to 'infinity'. Field lines will also open up in the presence of a powerful wind 
of particles emitted from the neutron star. Near the surface of the star, dipole magnetic field 
pressure is high enough to completely dominate the wind stresses. However, the magnetic 
energy density drops much faster with distance than that of a quasi-isotropic particle wind 
of kinetic luminosity Lp at infinity; thus, ignoring a transitional region between these two 
regimes, beyond a distance of the order ropcn ~ ''"o(-Bo'"o'^/2Lp)^/^ magnetic field lines will be 
'combed out' by the wind. Here, Bq is the value of the neutron star surface dipole magnetic 
field and Tq ~ 10 km, is the radius of the star (henceforth, we will use small r to denote 
spherical distances [from the center] , and capital it! for cyhndrical distances [from the axis] ) . 
The fraction of open field lines originates from a region of radial extent 

Rpc ~ Tq (ro/ropen) < Tq ( 1 ) 

around the axis, the so-called polar cap (this estimate is obtained for an undistorted dipole). 

An aligned magnetic rotator, the simplified geometry examined herein, spins down 
(even though an aligned magnetic rotator in vacuum does not), because a neutron star 
magnetosphere cannot be a true vacuum (Goldreich & Julian 1969). As is discussed in 
Contopoulos, Kazanas & Fendt (1999; hereafter CKF), a neutron star magnetosphere is 
spontaneously charged in order to support the charges and electric currents required in the 
reahstic solution. An electric current also flows, in a large scale electric circuit through the 
polar cap to infinity, closing along an equatorial current sheet discontinuity, which connects 
to the interface between open and closed field lines at the light cylinder, with the circuit 
closing across field lines along the polar cap. This electric current Ipc fiows along the 
magnetic field lines crossing the polar cap and generates the required spin-down magnetic 
torque on the neutron star. The neutron star spins down because the electromagnetic 
torque generated at the surface of the polar cap opposes its rotation. Inside the neutron 
star surface, this current flows horizontally towards the edge of the polar cap, where it 
flows out in a current sheet along the interface between open and closed held lines. The 
electric current flowing through the polar cap is (to a good approximation) distributed as 
/ oc \&(2 — "^/"^pc), where \& is the total amount of magnetic flux contained inside cylindrical 
radius R, and is the total amount of magnetic fiux which opens up to infinity. This 
is an exact expression for a magnetic (split) monopole, and CKF showed that it remains 
approximately valid even for a dipole. Since on the neutron star surface "if oc when 
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R< Rpc<t: ro, 
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liR) 




^ Rpc , 
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When this current flows horizontally in a layer of thickness h{R) across the polar cap, an 
electric current density 



J{R) - 



27rRh{R) 




JL 



(3) 



will flow horizontally, which, combined with the axial magnetic field -B^, threading the 
polar cap, generates an azimuthal Lorenz force per unit volume, f{R) = \J{R) x Bq. 
Integrating f{R) over the volume of the polar cap crust where the above electric current 
flows horizontally, and doubhng our result to account for the two (north/south) polar caps, 
we obtain the total electromagnetic torque 

2 



^JpcBoRpc- (4) 

We present two simple, physically equivalent, estimates of the electric current flowing in 
the magnetosphere. One is to consider particles (electrons/positrons) with Goldreich- Julian 
charge densities pcj ~ Bq/2'kRi^c flowing outwards at the speed of light from the polar cap. 
This naive estimate gives 

/,~.fl,P0.c= — (— j (5) 

Another equivalent, more physical, way to estimate the total amount of electric current 
flowing in the magnetosphere through the polar cap is to make the naive (and correct) 
assumption that, at the distance of the light cylinder the two magnetic fleld components 
(toroidal and poloidal) must be of the same order of magnitude, -B^|lc ~ Bp\i^c ■ This is 
indeed true in the force-free axisymmetric magnetosphere (without wind), since the light 
cylinder is the Alfven point (Li & Melrose 1994). In general, the Alfven point is some short 
distance inside the light cylinder. When the two fleld components are scaled back to the 
surface of the star at the edge of the polar cap, we obtain 

^ RlC \ f ropen 



Bp\pc = Bp\^c[^ ^Bo. (6) 

\ropcnJ \ ro J 

The structure of the fleld is dipole-like out to ropen and monopole-like out to the light 
cylinder. From Eqn (H) and the relation -B</,|lc ~ -BpIlc, we have 

^3 

B^lpc = Bo ° , and therefore (7) 

-'I'pc-'l'LC' open 
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B^l pcRpcC _ Bprpc / jtq \ f 

\ ' open / 

This is a very simple result, and shows that the two ways of looking at the problem are 
equivalent. 

Using the above relation for the polar cap current in equation we obtain the 
expression for the torque T associated with the above spin-down arguments. The 
corresponding energy loss rate due to the above torque is therefore E = T ■ Q, given more 
explicitly by 

where we have used Lp/incrl^^^ = i?(ropon)^/87r to obtain an expression for in 
terms of the wind luminosity, Lp, and dipole energy loss, E^,. Note that the standard 
dipole spin-down formula (modulo the different numerical factor in the denominator) is 
modified by the term (-RLc/''^open)^ which incorporates the effects of the "loading" of the 
magnetosphere with the outflowing wind. This expression, while of similar functional form, 
differs from that of TB98 because of a normalization error in this work (Thompson, priv. 
comm.), but agrees with a corrected expression given by Thompson et al. (1999). 

Using Eqn. (|) with values for P = Pigoe = 7.48 s, P = Asoe = 8.3 ■ IQ-^^ s s-\ i.e. 
those observed in SGR1806-20 and assuming the presence of a steady wind of luminosity 
L37 = Lp/10^'' erg s~^, we obtain an estimate of the surface magnetic field of the neutron 
star: 

Po - 3 X lO^^G (-^)" f^) L3//^ (10) 

where we have assumed ro = 10 km and / = lO^^gcm^. If one uses the observed value of 
the spin-down rate and the average value of the particle luminosity needed to account for 
the energetics of the nebula, then Bq is significantly below the estimate of Bq ~ 10^^ using 



Eqn (I) with Plc 



open- 



This modified spin-down law leads to exponential increase in the pulsar period instead 
of the power law increase associated with the purely dipole emission (this is easily seen from 
equating E'^, = —IQfl and integrating). One can thus estimate the age r of SGR1806-20 
through the relation 

P 



In 

2P 



LpP' 
47r2/P 



(11) 



which yields r ~ 11,800 yr for Lp = 10^^ erg s~^. Thus the steady wind model can naturally 
account for the fact that the age of the SNR GlO.0-0.3 is much larger than the characteristic 
dipole spin-down age, but with the penalty that the magnetar model must be abandoned. 
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3. COMBINED WIND AND DIPOLE SPIN-DOWN 

The expression for the magnetic field and characteristic age of the neutron star given 
above are only valid if the relativistic particle wind completely and continuously dominates 
the spin down of the star. If the particle wind fiow is either discontinuous or not dominant, 
then a more general description of the spin-down energy loss must be used. If the wind has 
instantaneous luminosity Lp during its times of activity and duty cycle Dp, defined as the 
fraction of total on-time, then the average energy loss from combined wind and dipole is 

E = -{mn) = Ed{1 - Dp) + E^Dp 

where we have used Eqn (^. The surface magnetic field may then be found as the solution 
to the quadratic equation, giving 



where. 



1/2- 



F(P,P) = ^1-(1+ \ ^ (14) 



l^vDl 



Bo ^ ^ ^ J (15) 



and E = A7r^I{P) /P'\ Note that when LpD^ < 4E(1 - Dp 

1/2 

27rV6(l - Dp)' 

which gives the standard dipole formula when Dp = 0. If LpDp ^ 4i?(l — Dp), Eqn (|T3|) 
gives the pure wind formula Eqn(|TU|) with Lp replaced by LpD^. 

We may also integrate Eqn (|l^) from the initial period Pq to the present period P 
to obtain the general expression for the neutron star characteristic age r. Assuming that 
Po « P, 

_ An^I ln[l-2/F(P,P)] 
- LpD^P^ F{P,P) ■ ^ ' 

This expression gives the usual characteristic age for dipole spin down, r = P/2P, in the 
hmit LpDp -C 4:E{1 — Dp) and Dp = 0. One must be careful in the limit LpDp 3> 4^(1 — Dp). 
If Dp is close to 1, then the first term in Eqn (|12]) should be dropped to give an expression 
for r which is the same as Eqn ([TT|) , again with Lp replaced by LpDp. 
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We now examine the consequences of the general expressions (0) and ([T^) for 
SGR1806-20, the SGR source for which we have the best estimate of the particle wind 
luminosity. Figure 1 shows the values of Bq and r computed from Eqns (JT^) and (JT^) for 
the measured P and P of SGR 1806-20 (Kouveliotou et al. 1998) as a function of the 
parameter LpD^/E, indicating the fractional wind contribution to the spin-down energy 
loss rate, assuming Dp <^1. For small LpD^/E, the curves approach their dipole radiation 
values of Bq ~ 10^^ G and r ~ 1500 yr. For LpD^/E ~ 0.1, Bq and r begin to depart from 
these values, with the magnetic field decreasing and the age increasing to connect smoothly 
to the wind-dominated solutions. We have seen in Section II that assuming continuous, 
wind dominated spin-down can give a characteristic age that agrees with the age (r ~ 10"^ 
yr) of the plerion surrounding SGR1806-20, but that the derived magnetic field drops into 
the range of normal radio pulsars. In such a case, however, the free energy associated with 
the magnetic field decay is not sufficient to account for the observed luminosity, and an 
alternative free energy source must be considered. 

The parameter LpD^ may be estimated for SGR 1806-20 from its observed burst 
characteristics. In general, we can write the particle luminosity associated with a burst 
as Lp = E^e^^ '^'^w^^ where E^ is the 7-ray burst energy, is the conversion efficiency of 
particle energy to 7-rays and Ar^ is the duration of the wind outflow following the burst. 
If T is the average time between bursts, then the wind duty cycle is Dp = At^/T, and 

LpDl = E,t:^^Ar^T-\ (17) 

In addition, the requirement that the X-ray nebula (Murakami et al. 1994) is 
powered by the aggregate of the bursts' wind outflows leads to the condition 
LpDp = 10^^erg£'4o(10"^/e^)/T = lO^^ergs"^ (lO"^///), where r] is the conversion 
efficiency of particle luminosity to nebular emission and E^q = E^/ 10'^^ erg. 
For the multiple SGR bursts from SGR 1806-20, Eqn (|l^) and the above 
requirement on LpDp, we find i?4o (10~^/e^) = (T^g/j/IO^ s)(10^^/?7), giving 
LpDl = lO^^ergs-i (Ar^/102s)(10^s/T5GR)(10"7r7). Likewise, for the giant bursts, 
we have LpD^ = 3 x lO^^ergs^^ (Ar^/10^ s)(10 yr/rG)(10-V^)- From Eqn (|13D, Eqn 



and Figure 1, the conflicting goals of preserving the magnetar model (i.e. Bq ~ 10^^ G) and 
bringing the characteristic age within a factor of 2 of the 10^ yr age of GlO.0-0.3, may be 
satisfied with LpD^/E ~ 10 — 100. Since E = 8 x lO^^ergs^^, the duration of the particle 
outflow must be much larger than the 7-ray burst duration and the wind flow duty cycle 



must be. Dp ~ 0.008 - 0.08. 
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4. DISCUSSION 

The results of our analysis leave two alternatives for interpreting the spin down of 
SGRs, given the present limited data. The first assumes a continuous wind outfiow at a 
luminosity sufficient to yield a characteristic age in agreement with that of the surrounding 
SNR. We have shown that in the case of SGR1806-20, this assumption leads to a surface 
magnetic field of only Bq — Sx 10^^ G, well below the magnetar range (> 10^^ G). However, 
this alternative requires a source of free energy other that the field decay to power both 
the nebular emission and the SGR bursts. The second option assumes an episodic (or at 
least variable) wind outflow such that the average wind luminosity is sufficient to power 
the nebular emission. This allows a range of combinations of surface field and age and 
depends on the wind duty cycle. We show that, in the case of SGR1806-20, it is possible 
to accommodate both a characteristic age r ~ 7500 yr, consistent with the estimated SNR 
age of ~ 10^ yr, and a magnetar model {Bq — lO^^G). One should then observe a sudden 
increase in the period derivative following SGR bursts. Evidence for such an increase was 
seen following the bursting activity of June - August 1998 from SGR1900+14 (Marsden et 
al. 1999). These options assume that gravitational radiation did not play a significant role 
in early spin-down evolution of the star, but it is unlikely to have made a large difference in 
the characteristic age. 

There are a number of arguments in favor of the magnetar model for SGRs, most of 
which have been discussed by Thompson & Duncan (1995) and Baring & Harding (1998). 
An additional argument is that the pure dipole fields of AXPs, which are spinning down 
smoothly (Gotthelf et al. 1999) and have much lower luminosity wind flows (if any), lie in 
the magnetar range. We suggest that detailed monitoring of the spin periods of the SGRs, 
to search for variations in the period derivative, can measure or place limits on the duty 
cycle of particle outflows and thus determine whether a magnetar model for these sources 
is viable. 

We thank Rob Duncan, Chris Thompson, Peter Woods and the referee Cole Miller for 
valuable commments and discussions. 
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Fig. 1. — Surface dipole magnetic field Bq and characteristic age r as a function of the 
parameter LpD^/E, from the general expressions Eqns (0) and ([T6| ) derived from the 
combined dipole and wind spin-down formula, Eqn (0). 



